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We constrain the Higgs boson (Yukawa) coupling to quarks in the first two generations in the 
H ^ ZZ final states. Deviation of these couplings from the Standard Model values leads to change 
in the Higgs boson width and in the cross sections of relevant processes. In the Higgs boson resonance 
region, an increased light Yukawa coupling leads to an increased Higgs boson width, which in turn 
leads to a decreased cross section. In the off-shell region, increased Yukawa couplings result in an 
enhancement of the Higgs boson signal through qq annihilation. With the assumption of scaling 
one Yukawa coupling at a time, this study is conceptually simple and yields results with the same 
order of magnitude as the tightest in the literature. The study is based on results published by 
the CMS experiment at the LHC in 2014, corresponding to integrated luminosities of 5.1 fb~^ at a 
centre-of-mass energy y/s = 7TeV and 19.7fb“^ at 8TeV. 


Since the discovery of a Higgs boson with a mass 
of around 125.6 GeV at the Large Hadron Collider 
(LHC) [I|, B] , measurements of its properties have shown 
consistency with the Standard Model (SM) expectations 
within the uncertainties il. Assuming SM, the gluon 
fusion via closed quark loop dominates the Higgs bo¬ 
son production, because of large gluon-gluon luminos¬ 
ity and large mass of the top quark. Also due to rel¬ 
atively large mass of the b quark, the Higgs boson de¬ 
cays into a pair of b quarks most of the time. While 
experimental analyses have been performed on the inter¬ 
actions between the Hi ggs boson and heavy quarks [i- 
d, as well as leptons [l3l - [l7l |- no experimental results 
have been presented on the Higgs Yukawa coupling to 
the light, namely, u, d, and s quarks. This is not surpris¬ 
ing, because in SM, (1) the small masses of the u and d 
quarks make their Yukawa couplings to the Higgs boson 
weak, with the branching fraction (B) of the Higgs bo¬ 
son decaying to it or d quark pair being < 10“®; and 
(2) while B{H —>■ ss) ^ 2.4 x 10“"^ (comparable to 
B{H —>• /!+/!“)) and B{H —>• cc) ^ 2.9 x 10“^ (com¬ 
parable to B{H ^ ZZ)), these decaying processes are 
difficult to observe without efficient quark flavor tagging. 
It is worth noting that in Ref. [l^, the upper limits of 
B{H —>• and B[H —>• e+e“) have been set to 

be 0.0016 and 0.0019 respectively, where the latter is 
Ri 3.7 X 10® times the SM value. It is also worth not¬ 
ing that phenomenological studies do exist on constrain¬ 
ing light Yukawa couplings of light quarks. For example, 
in Ref. [l^, depending on the analysis performed, the 
upper limit of the Yukawa coupling between the Higgs 
boson and the c quark can be set as low as 5, 6.2 times 
the SM value. Another example is Ref. [2, in which 
Higgs-boson-mediated production of vector meson in as¬ 
sociation with a vector boson is used to constrain Yukawa 
couplings of u, d, and s quarks. Via a global fit and de¬ 
pending on how the couplings are allowed to vary, the 
upper limit of these couplings are found to be close to 
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the SM Yukawa coupling of the b quark, whose numeri¬ 
cal indications will be summarized in a later table. 

In this study we attempt to constrain the Yukawa cou¬ 
pling between the Higgs boson and quarks in the first two 
generations. The Yukawa Lagrangian after electroweak 
symmetric breaking is 

^Yukawa = - X - ~ff^^ (1) 

t; 

where the summation is over fermion flavors. Relaxing 
the coupling constants, the deviation from SM considered 
in this study is written as 

A/lvukawa = - 1 )^//H, ( 2 ) 

where f = u,d, s, c, and the scaling factors c/ are real and 
can take both positive and negative values. The masses 
of the quarks are set according to the 2014 Particle Data 
Group summery table . As we keep the couplings the 
Higgs boson to other particles SM, the deviation of c/ 
from 1 leads to change in the Higgs boson width, and in 
the cross sections of processes involving Yukawa interac¬ 
tions. 

Combining the direct measurement with 77 and Af fi¬ 
nal states, the CMS experiment has set an upper limit 
for the Higgs boson width at 1.7 GeV at a 95% conh- 
dence level (CL) [^. This can be translated to an upper 
limit for each |c/| by adding its contribution to the Higgs 
boson width predicted by SM. In Table [H we list the up¬ 
per limits on |c/| due to this argument, the calculation is 
performed at leading order (LO), with the Higgs boson 
mass niH = 125.6 GeV and the corresponding SM width 
r™ = 4.15 MeV. As reference, we list the upper limits 
by requiring the theory being perturbative, namely, 

< 0(1). (3) 

In addition, we list the upper limits by Ref. [l^ and 
Ref. [1^, as some of the best constraints placed thus far. 
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\Cu 1 

\cd\ 

|Cs| 

|Cc| 

Perturbation 

< 1.1 X 10® 

< 5.1 X 10“^ 

< 2600 

< 190 

Th < 1.7 GeV 

< 4.9 X 10^ 

< 2.4 X 10“^ 

< 1200 

^ 88 

Ref. [1^ 

Ref. [IJ 

2100 - 2800 

930 - 1400 

35-70 

^ 6.2 


TABLE I: 95% CL upper limits of scaling factors |c/|, 
due to Higgs boson width direct measurement; in 
comparison with those by requiriM the theory 
perturbative, and those by Ref. [3 and Ref. [l^. In 
the Standard Model, c/ = 1. 


In the remainder of this study we explore constraints 
on c/ from the production of the Higgs boson, which de¬ 
cays into a pair of Z bosons. In the standard model, 
the production of the Higgs boson is dominated by gluon 
fusion with t and b loops (Fig. [TJa)). The dominant con¬ 
tinuum background is the quark-initiated ZZ production 
(Fig-HKb)), accompanied by gluon-initiated ZZ produc¬ 
tion (Fig. [TJc)). The subdominant production mecha¬ 
nism of the Higgs boson is vector boson fusion (VBF, 
Fig. md)), which contributes about 7% to the Higgs bo¬ 
son production in the resonance region, and about 10% 
in the mzz > ‘i-'mz region. 
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FIG. 1: LO contributions to the main ZZ production 
processes: (a) Higgs-mediated gg production; (b) 
quark-initiated background production; (c) 
gluon-initiated background production; and (d) 
Higgs-mediated VBF production. 


As the Yukawa couplings change with c/, additional 
contributions from the Higgs-mediated quark annihila¬ 
tion (Fig.[2Ka, b)) and gluon fusion with light quark loops 
(Fig. (He)) are taken into account in this study. While we 
understand that a large |c/| could make a difference in 
the VBF type diagram, by having Higgs boson in place of 
the weak bosons (Fig. EJd)), we neglect its contribution 
in this study because of its distinct kinematic characteris¬ 
tics, particularly the angular correlation between the two 
jets, which will allow suppression (See, e.g.. Ref. [^L 
For a typical process with the Higgs boson created in 
the s-channel, the cross section in the resonance and off- 



FIG. 2: Additional LO contributions to the ZZ 
production due to b and light quarks: (a) 
Higgs-mediated bb annihilation; (b) Higgs-mediated 
light-quark annihilation; (c) Higgs-mediated gg 
production; and (d) H ZZ production via Higgs 
boson fusion. 


shell region behaves as 


resonance 


9ah 9hb 

Th 


and 


off-shell 


2 2 
9AH 9hb^ 


(4) 

respectively, where gAH (gHB) is the Higgs boson cou¬ 
pling to the initial (final) state. The resonance and off- 
shell regions offer two distinct sources of information on 
the Yukawa couplings. In the resonance region, the gluon 
fusion is always the dominant mechanism of Higgs bo¬ 
son production, even at large scaling factors |c/|, as a 
result of the large gluon-gluon luminosity. While each 
scaled Yukawa coupling makes its contribution in the 
gg ^ H closed quark loop, its contribution to the Higgs 
boson width affects the production overwhelmingly at 
large |c/|, reducing the cross section to near zero. This 
feature can be used to constrain c/ by requiring |c/| being 
small enough to allow consistency with experimental ob¬ 
servations. In the off-shell region, the production cross 
section of the Higgs boson, which decays into ZZ, re¬ 
ceives enhancement at mzz ^ “^mz, where the invariant 
mass of the Hig gs boson allows both Z bosons become 
on-shell [2^, |24| In addition, the parton luminosities of 
qq, particularly uu and dd, are less dominated by that of 
gg (See, e.g., Ref. [l^). As a result, the production of the 
Higgs boson becomes dominated by the qq annihilation 
(Fig. [IKb)) at large |c/|, and the cross section increases 
with |c/p. Therefore the measured off-shell cross section 
may be used to further constrain c/. 


We now calculate the cross section oi ZZ production 
with four-lepton (4£, i = e, /i) final states in proton- 
proton collisions at centre-of-mass energy y/s = 7TeV 
and 8TeV, as a function of each individual Cf, and com¬ 
pare to the signal strength g-ggH reported by CMS in 
Ref. @,[21, which is based on integrated luminosities of 
5.1 fb”^ at Vs = 7TeV and 19.7fb"^ at 8TeV. The 
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signal strength is defined by the relation, 

^obs. _ ,, SM , /-- SM , SM 

^gg^ZZ t^ggH ^signal yt^ggH ^intf. ' ^bkg.5 

where <y°f^zz^ and a™ are the to¬ 

tal gluon-initiated ZZ cross section observed, Standard 
Model predictions for the Higgs boson signal, signal- 
background interference, and background, respectively. 
While our analysis largely involves Higgs boson pro¬ 
duction by qq annihilation, its indistinguishability from 
gluon fusion allows us to base the analysis on fJ^ggH- In 
each calculation, one c/ is varied in the range set by 
Table |T] while others are kept at 1 (SM value). The reso¬ 
nance region is defined as 120.5 GeV < 7714 ^ < 130.6 GeV 
in our calculation, which does not necessarily agree with 
the CMS dehnition. The off-shell region is defined as 
220 GeV < 7714 ^ < 800 GeV, as adopted by CMS in 
Ref. [Hi. 

The gluon-initiated processes are calculated as follows. 
The contribution from Fig. [IJa),[21(c), and their interfer¬ 
ence with Fig. (iKc) is calculated using MCFM 6.8 (H 
with NPROC = 128 — 130 at loop-induced leading order 
in perturbative quantum chromodynamics (QCD). The 
contribution from Fig. He) is implemented by adding 
codes that are parallel to those that calculate the t 
and b loops. The contribution from Fig. [T] (c) is cal¬ 
culated with NPROC = 81 at loop-induced leading order. 
For simplicity, the quark mass evolution is accounted as 
part of the uncertainty in cross section. The cross sec¬ 
tions of the gluon-initiated processes are scaled by the 
same m^^-dependent correction factors to the LO cross 
section {K factors) applied in Ref. [Hi with next-to- 
next-to-leading order and next-to-next-to-lea ding loga¬ 
rithms accuracy for the total cross section [28l - [30l| - The 
QCD renormalization and factorization scales are set to 
Mr = M/ = (dynamic scales) and MSTW2008 

parton distribution functions (PDFs) [H are used. 

The quark-initiated processes are calculated using 
MadGraphS V2.2.3 [H at LO, with a modified SM 
model that includes the Yukawa couplings for all the 
quark species. Interferences between Fig. Hb) and 
Fig. IJJa, b) are given full treatment. For simplicity, the 
quark mass evolution is accounted by scaling down the 
cross sections by a factor of one half. The cross sec¬ 
tions calculated at LO are scaled by a K factor, which is 
the ratio between the NLO cross section of Fig. [TKb) ob¬ 
tained by MCFM and the LO one by MadGraphS. The 
cross section is further scaled hy a K factor due to NLO 
electroweak corrections [H IH i which predicts negative 
and m^^-dependent corrections to the qq ^ ZZ pro¬ 
cess for on-shell Z boson pairs, fir = = mzzj^ and 

MSTW2008 PDFs are used. 

We employ the CMS selection cuts [Hi requiring 
Pa.m > 5 GeV, > 7 GeV, |r;^| < 2.4, < 2.5, 

'mi+e- > 4 GeV, M 41 > 100 GeV. In addition, the trans¬ 
verse momentum of the hardest (next-to-hardest) lepton 
should be larger than 20 (10) GeV, the invariant mass of 
a pair of same-flavor leptons closest to the Z boson mass 


should be in the interval 40 GeV < m^+(- <120 GeV 
and the invariant mass of the other pair should be in the 
interval 12 GeV < m^+i- < 120 GeV. 

In the resonance region, the cross section of qq anni¬ 
hilation is dependent on c/, but its contribution to the 
total cross section is tiny due to low parton luminosities 
of qq. For example, even at c„ = 5 x 10^, this cross sec¬ 
tion is only at a few percent level of that of the SM gluon 
fusion. The strongest dependence on c/ is with the gluon 
fusion process. Fig. [TJa) and Fig.j^Kc), combined with its 
interference with the continuum background. This cross 
section decreases as |c„| increases, because of its 1 /Th 
dependence, as suggested in Eq. O- At sufficiently large 
|cu|, the signal is washed out, and the gluon-initiated pro¬ 
cesses reduce to the continuum background production of 
ZZ. 

In Fig. |3l we combine the results in the resonance re¬ 
gion, and present the expected signal strength, translated 
from cross section by Eq. ®i as function of each indi¬ 
vidual c/, and compare to the signal strength for gluon 
fusion UggH = 0.85to;i6, reported by the CMS experi¬ 
ment In Table m we summarize the 95% CL limits 
for each c/. We notice here a slight asymmetry about 
Cf = 0 due to the sign of the interference in the gluon- 
initiated processes. 

-3300 < Cu < 3200 

-2000 < Cd < 1900 

-130 < Cs < 125 

-11 ^ Cc < 9 

TABLE H: 95% CL limits of scaling factors c/, by the 
observation of production of the Higgs boson in the 
resonance and its decay to the ZZ final states. 

In the off-shell region, the Standard Model expects the 
sum of Higgs boson signal and its interference with the 
continuum background to be slightly negative. While this 
sum is dependent on |c/|, the total cross section becomes 
dominated by the quark-initiated process at large |c/|. 
Although the sensitivity of the off-shell cross section is 
not as high as that of the gluon-initiated process in the 
resonance region for relatively small |c/|, sufficiently large 
|c/| give rise to a departure from the number of events 
observed by CMS. 

In Eig. m we compare the expected off-shell signal 
strength as functions of c/ with the one estimated from 
the result published by CMS in Ref. [H IH ■ 
here that background subtracted number of events, in¬ 
stead of signal strength, are presented in Fig. H) The 
number of events are translated from signal strength by 
Eq. (I5|). Moreover, since the observation of the Higgs bo¬ 
son signal in the resonance is a well established fact, the 
analysis respects this fact by scaling the off-shell signal 
strength in such a way that the resonance signal strength 
is fixed. In other word, from the analysis in the resonance 
region, we learned that signal strength decreases with in¬ 
creasing |c/|, and therefore in the off-shell analysis, the 
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FIG. 3: Number of ZZ —>■ 4£ events expected in the resonance region (105.6 GeV < 7714 ^ < 140.6 GeV) as functions 
of Cf (Icr and 2a uncertainties in green and yellow bands), in comparison with number of events observed by GMS 
( 1(7 and 2a uncertainties in red and pink bands), with 5.1 fb~^ proton-proton collisions at yfs = 7TeV and 19.7fb~^ 
at 8 TeV. See text for details of calculation. 


signal strength is scaled accordingly so that for any |c/| 
the resonance signal strength remains 0.85. In Table nni 
we summarize the 95% GL limits for each c/ set by the 
analysis in the off-shell region. 

\cu\ < 1.5 X 10^ 

|cd| < 8600 

|cs| < 1100 

|cc| < no 

TABLE III: 95% CL upper limits of scaling factors |c/|, 
by the observation of off-shell production of the Higgs 
boson and its decay to the ZZ final states. 

From the analysis in the off-shell region, we have ob¬ 
tained upper limits for \cu\ and \cd\ that are over twice 
tighter than those due to the Higgs boson width direct 
measurement, and slightly tighter results for |cs| and |cc|. 
The better performance of this analysis on \cu\ and \cd\_ 
is due to the higher parton luminosities of uu and dd 
than those of ss and cc. Furthermore, at a hxed energy 
of proton-proton collision, although parton luminosities 


decrease in general as the center-of-mass energy of the 
colliding partons increases, the rates of decreasing for uu 
and dd luminosities are slower than those for ss and cc, 
which are still slower than that of gg (See, e.g.. Ref. [^1. 
This suggests an improvement of this analysis as we ex¬ 
plore the higher invariant mass region. As an illustra¬ 
tion, we perform the analysis with 7714 ^ > 1200 GeV at 
High-Luminosity LHC (HL-LHC), where an integrated 
luminosity of 3000 fb~^ is delivered at ^/s = 14 TeV. As 
shown in Fig. O a factor of ^ 3 improvement over the 
current limits can be expected for \cu\ and \cd\, as well 
as a factor of ~ 2 for |cs| and |cc|. 

The analysis in the off-shell region may be further im¬ 
proved in two ways. The first is to improve the statistics 
and include the WW —>■ 2£2iy final states. The second 
is to suppress the ZZ {WW) continuum background by 
emploYing a matrix-element-based method, as done in 
Ref. 1^ . I26lp1. We also note that, due to the lack of ex- 


^ See, e.g. Ref. l22l l36l for more details on the signal-background 
separation for H ^ ZZ and H —7- WW processes using the 
matrix element likelihood analysis. 
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FIG. 4: Number of ZZ —>■ 4£ events due to Higgs boson (signal and interference) expected in the off-shell 
(220 GeV < m 4 £ < 800 GeV) region as functions of c/ (Icr and 2cr uncertainties in green and yellow bands), in 
comparison with number of events observed by CMS (Icr and 2cr uncertainties in red and pink bands), with 5.1 fb~ 
proton-proton collisions at = 7TeV and 19.7fb“^ at 8TeV. See text for details. 


perimental access, we were unable to perform the analy¬ 
sis with combined CMS and ATLAS results. We believe, 
however, an analysis based on ATLAS result would yield 
very similar constraints; and for the purpose of demon¬ 
strating our analysis methods, it suffices to use CMS re¬ 
sults alone. 

In conclusion, we suggested that the Yukawa coupling 
between the Higgs boson and light quarks can be con¬ 
strained by comparing the signal strength of Higgs boson 
production as functions of scaling factors of the Yukawa 
couplings with the value measured by the LHC experi¬ 
ments. The tightest constraints are set in the resonance 
region, which are listed in Table |TT1 With the assump¬ 
tion of scaling one Yukawa coupling at a time, these con¬ 
straints are at the same order of magnitude as the best 
phenomenological results in the literature [Ulil. While 
the analysis performed in the Higgs off-shell region does 
not place as tight limits on the scaling factor, it places 


independent constraints. In addition, we have demon¬ 
strated that the analysis performed in a higher invari¬ 
ant mass region can receive improvement by taking ad¬ 
vantage of the behavior of parton luminosity functions. 
While we believe that our analysis is sufficiently accu¬ 
rate for order-of-magnitude estimates, the present study 
is crude and ignores the many theoretical as well as exper¬ 
imental details. Therefore, it will be best if experimental 
collaborations perform a detailed analysis. 
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FIG. 5: Number of ZZ —>■ 4£ events due to Higgs boson (signal and interference) expected in the off-shell 
(to 4 ^ > 1200 GeV) region as functions of c/ (Icr and 2cr uncertainties in green and yellow bands), in comparison 
with number of events expected by HL-LHC (Icr and 2cr uncertainties in red and pink bands), with 3000fb”^ 
proton-proton collisions at -Js = 14 TeV. See text for details. 
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